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BENDING AND STRETCHING UNIT VECTOR 
FIELDS IN EUCLIDEAN AND HYPERBOLIC 3-SPACE 



C. M. Wood 

For Esther. 



• Abstract. New examples of harmonic unit vector fields on hyperbolic 3-space are con- 

structed by exploiting the reduction of symmetry arising from the foliation by horospheres. 
This is compared and contrasted with the analogous construction in Euclidean 3-space, 
using a foliation by planes, which produces some new examples of harmonic maps from 
3-dimensional Euclidean domains to the 2-sphere. Finally, the harmonic unit vector field 
tangent to a parallel family of hyperbolic geodesies is shown to be unstable, by constructing 
a class of compactly supported energy-decreasing variations. All examples considered have 
infinite total bending. 
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\q . 1. Introduction 

O 

A smooth unit vector field a on an orientable Riemannian manifold (M, g) is said to 

a ■ 

be harmonic if a is a harmonic section of the unit tangent bundle of M [17, 18]: 

> 

V*Va = |Va| 2 a, (1.1) 

a: 

where V is the Levi-Civita connection and V*V = — TrV 2 . The coefficient |Vcr| 2 , 
which is (twice) the vertical energy density of a, is sometimes also called the bending 
of a [17], and (1.1) are the Euler-Lagrange equations for the associated functional on 
unit fields. The study of harmonic unit fields has become popular in recent years, as 
witnessed by the bibliography of [8]. There are a number of results specific to the 3- 
dimensional case (see [10, 11, 12, 14] amongst others): for example, if a is a unit Killing 
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2 C. M. WOOD 

vector, or the characteristic vector of a contact metric structure, then a is harmonic if and 
only if a is (pointwise) an eigenvector of the Ricci operator. The same characterization 
applies when a is tangent to the fibres of a harmonic morphism from M 3 to a surface 
[2]. Furthermore, in space forms of all dimensions, radial unit vector fields about a 
point, or a totally geodesically submanifold, are harmonic [3]. These lead to a number 
of "natural" examples, such as a classification of left-invariant harmonic unit fields on 
3-dimensional unimodular Lie groups, and consequent classification (including stability) 
of invariant harmonic unit fields on compact quotients. A unified account of the local 
and global classification of harmonic morphisms from 3-dimensional space forms can be 
found in [2] . In this paper a number of "less natural" examples of harmonic unit fields 
are constructed, when M is an open subset of M 3 or H 3 . In the Euclidean case we may 
view a as a map into S 2 , and (1.1) is then equivalent to the harmonic map equations 
[7]; thus, we are constructing harmonic maps M — > S 2 . However since this is no longer 
the case when M C H 3 we will work with equations (1.1) throughout. The basic idea is 
to simplify (1.1) by breaking the symmetry of M and asking a to respect this additional 
structure. The symmetry reduction is guided by the availability of a reasonable supply 
of simple examples of harmonic unit fields which satisfy the corresponding constraint. 
In the Euclidean case we simply ask a to be planar; that is, tangent to a foliation by 
parallel planes. We may then (locally) factorise a: 

where w(t) = e lt and S 1 is included in S 2 as the equator cut out by the foliation, and this 
reduces (1.1) to the linear equation Au = 0. It is then quite straightforward to restore 
the full symmetry in order to construct non-constrained harmonic a (Proposition 1). One 
equation in the resultant system of nonlinear PDEs may be viewed as a generalization of 
the elliptic 3-dimensional sine-Gordon equation. In §3 we use this approach to construct 
a family of continuous maps 1R 3 — > S 2 , parametrized by an open Mobius band, all of 
which are harmonic (hence smooth) on the complement of a line, and whose energy 
density extends to a bounded continuous function on M. 3 . 

In §4 we consider unit vector fields in H 3 which are horospherical; that is, tangent to a 
foliation of M by horospheres. Working with the upper half space model, equation (1.1) 
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then reduces to a linear second order equation (viz. the hyperbolic Laplace equation) with 
a nonlinear first order constraint (Proposition 2). The latter cuts down the availability 
of solutions quite dramatically compared with the Euclidean case. Nevertheless, we are 
able to construct a smooth 1-parameter family of harmonic unit fields on H 3 , no member 
of which is invariant under any discrete subgroup of isometries with compact quotient, 
but whose energies are arbitrarily close to the minimum for horospherical unit fields, this 
being achieved only by the invariant unit fields. In §5 we examine the stability of unit 
fields a whose integral curves are geodesies which are parallel in the hyperbolic sense ie. 
convergent to a common point on the ideal boundary of H 3 . We refer to such fields as 
H-parallel. They are harmonic [9], and their variation fields are precisely the compactly 
supported horospherical vector fields, for the orthogonal horosphere foliation. In contrast 
to the Hopf vector field on the 3-dimensional sphere [4, 14], we show that a is unstable 
in the following way: there exists R u > such that for any geodesic ball B R of radius 
R > R u it is possible to construct a compactly supported energy-decreasing variation of 
a whose support lies inside B R . We summarise this by saying that a is stretchable in all 
sufficiently large geodesic balls. 

It is a pleasure to thank Evgeny Sklyanin, and the referee, for their perceptive and 
helpful comments. 

2. Generalities, and Euclidean Preliminaries 

We first record some general identities in Riemannian geometry related to equations 
(1.1). If X is a vector field on M and /: M — > R is a smooth function, where M is any 
Riemannian manifold, then: 



V*V(/X) = /V*VX - 2X7 vf X + AfX, 



(2.1) 



where V/ is the gradient vector, and Af is the Laplacian: 



Af = -TrVdf = -div(V/). 



Note our choice of sign for Af. Furthermore if g: R — > R is smooth then: 



V(<7(/)) = g'(f)Vf, 

= g\f)Af - g"(f)\Vf 



(2.2) 



(2.3) 
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Suppose now that M is a domain of Euclidean IR 3 , and a is a smooth unit vector field 
on M. Let (£1,^2, £3) be a right-handed orthonormal basis of R 3 , and let L denote the 
line generated by £3, and P the plane spanned by £1 and £ 2 . Let T be the foliation of 
M by planes parallel to P, and let J be the complex structure in T compatible with the 
induced orientation. Thus £2 = </£i- If cr(x) 7^ ±£3 for all x E M then a may be written 
in polar form: 

o = cos -U sin "u £x + sinwsini^ + cosf £3, (2.4) 

where u, v are smooth IR- valued functions. In general, u is only locally defined, any two 
differing by an integer multiple of 2-k. The equatorial part £ of a is then the unit field: 

£ = cosu^i + sinw^. 

Imposing the harmonic section equations (1.1) on a yields an over deter mined, but con- 
sistent, system of nonlinear PDEs for u and v. 

Proposition 1. Let a be a smooth unit vector field on a domain of Euclidean 3-space, 
expressible in polar form (2.4). Then the bending of a is: 

|Va| 2 = \Vv\ 2 + sin 2 w|Vw| 2 , (2.5) 

and a is harmonic if and only if: 

Au = 2(Vu»Vv)cotv, (2.6) 
Av = — I Vu\ 2 cos w sin (2.7) 

Proof. In terms of its equatorial part £, we have: 

a = sini;£ + cos?; £3. 

Now: 

V x £ = — (X.u) sin?/£i + (X.u) cos -u£ 2 = (X.u)J£, 

and therefore: 



V x cr = (X.v)cosv^ + (X.u) sinv J£ — (X.v) sinv £3, 
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V*V£ = A(cosw)£i + A(sinu)£ 2 

= (— sin it Am + cos u | V-u| 2 )£i + (cos u Aw + sin-u | Vw| 2 )^2 
= |V-u| 2 £ + Au Jf, 

and therefore: 

V*Va = sinvV*V^ - 2cosuV Vi; f + A(sinv)£ + A(cosi>)£ 3 

= sinv(|Vw| 2 ^ + AuJf) - 2(Vu«Vv)cosu J£ 

+ (cosvAw + sinw|Vw| 2 )^ + (— sinwAv + cosd | Vi>| 2 )£3 

= |Vcr|V+ (Av + | Vu\ 2 casvsm.v)(oosv£, — sini;£ 3 ) 
+ (Au sinv — 2(Vu»Vv) cos v) J£, 

where the final equation is obtained using (2.5). The result follows on comparison with 
(1.1), since (£, J£, £ 3 ) is an orthonormal frame. □ 

Remarks. 

(1) If |Vu| 2 is constant, say then (2.7) is the 3-dimensional elliptic sine- Gordon 
equation [5]: 

AV = —ji sin V, where V = 2v. 

(2) We say that a is planar if it is tangent to T . In this case v = 7r/2, the bending of 
a is | Vw| 2 , and cr is harmonic if and only if Aw = 0. 

(3) The natural C-action on J 7 , arising from its complex structure, induces a circle- 
action on planar unit fields. This extends to a circle-action on unit fields in polar form, 
when applied to their equatorial part; thus, the action of e lt on a is the unit field at with 
u t = u + 1. It follows from Proposition 1 that if a is harmonic then a t is a 1-parameter 
variation through harmonic unit fields, all of which have identical bending, but which in 
general are non-congruent. 

We illustrate Proposition 1 by first showing how it may be used to reconstruct some 
known examples of harmonic unit fields [3] . These will provide useful points of comparison 
for our forthcoming new examples. 
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Example 1. Let (r, 9, z) be cylindrical polar coordinates with respect to (£i, & , £3), and 
let u = 9 and v = tv/2. Then Au = 0, and the corresponding harmonic unit field is a = r , 
the unit radial field about L. The 1-parameter variation a t of a includes 0^/2 = 0, the 
unit vector in the ^-direction. All the a t are planar, with domain M = IR 3 \ L, and 
unbounded bending: 



r.- t | 2 = |w| 2 = \e/r\ 2 = 1. 



Example 2. Let (R,9,(f>) denote spherical polar coordinates relative to (£1,^27 £3), and 
let u = 9 and v = (p. Then Au = 0, and: 

Vw = 0/r, Vv = 4>/R. 

From the Laplacian in spherical polar coordinates: 

119, 1 sine/) cos 4> 

Av = — --7T— — -— -smffl = ——7 cot ( 



R 2 sin <j) d(j) R 2 r 2 

Therefore both (2.6) and (2.7) hold. The corresponding harmonic unit field is a = R, the 
unit radial field about the origin. The equatorial part £t of at is precisely the 1-parameter 
planar family constructed in Example 1; however, a t now extends to a harmonic unit 
field on R s \ {0}. The a t again have unbounded bending: 

In this case, a is tangent to the fibres of a harmonic morphism; indeed, when regarded as 
a map io S 2 , a is a. harmonic morphism. The existence of a variation through harmonic 
maps of identical energy density, none of which is a morphism, is an interesting feature 
of this example. 

3. New examples in Euclidean space 

Let £ denote any member of the 1-parameter family of planar harmonic unit fields 
from Example 1, defined on M = R 3 \ L. By a harmonic desingularization of £ we mean 
a harmonic unit field a on a strictly larger domain M satisfying: 



(3.1) 
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where A: M — > R is a smooth function. Example 2 may be regarded as a harmonic 
desingularization of £ over M = R 3 \ {0}. We now attempt to construct a harmonic 
desingularization on M = R 3 . Continuity of the left hand side of (3.1) implies that A 
vanishes on L. Now the derivative of A£ is: 

d\ ® i + \d9® J£, 

and since the derivative of the left hand side of (3.1) is continuous it follows that dX also 
vanishes on L. Successive differentiations of A£ lead to the conclusion that all derivatives 
of A vanish on L. Now, by the regularity of harmonic maps [6], a is real analytic when 
viewed as a function R 3 — > 5" 2 . Hence the left hand side of (3.1) is also real analytic, 
viewed as a function R 3 — > R 2 , and since all its /c-jets vanish on L, it is constant. Because 
a = £3 on L it follows that a = £3 identically, and the only harmonic desingularization 
of £ is therefore trivial: A = 0. 

We now show how to construct examples of continuous unit fields a on R 3 which 
satisfy (3.1) nontrivially and are harmonic on M = R 3 \ L. Expressing a in polar form 
(2.4) on M, it follows from (2.6) that 6»Vv = 0. Thus v is rotationally symmetric 
about L, and we may write v = v(r, z) in terms of cylindrical polar coordinates (r, 6, z). 
Suppose for simplicity that v = v(r). Then: 

-Av = v" + — , 
r 

and (2.7) reduces to the following nonlinear ODE with a singularity at r = 0: 

r 2 V" + rV = sinF, where V = 2v. (3.2) 
In logarithmic coordinates t = log r, (3.2) transforms to: 

V" = sinV, 

and putting x = tv — V yields the classical pendulum equation: 

x" + sinx = 0. (3.3) 

Continuity of a along L requires v(0+) = 0, which translates to x(— 00) = it. Furthermore 
if y = x' then: 

y{— 00) = — lim V'(s) = limrV'(r) = 0. 
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So in the familiar phase portrait of (3.3) (see for example [l,pp. 90-92]) the relevant 
trajectories are those on the separatrix which exit the critical point (tt, 0). These trajec- 
tories approach the neighbouring critical points (— tt, 0) and (37r, 0) as t — > oo, and we 
deduce that v = V/2 increases (resp. decreases) monotonically to tt (resp. — tt) as r — > oo. 
Thus as r — > oo the unit field cr approaches —£3, whereas a = £3 on L. The equation of 
the relevant piece of the separatrix in the phase plane is: 

y = — 2cos(x/2), — tt < x < 3tt. 

This first order ODE for x has solution: 

{c(sec(;r/2) + tan(x/2)), — tt < x < tt, 
— c(sec(x/2) + tan(x/2)), tt < x < Stt, 
where c > is a constant of integration, and it follows that: 

v'(0+) = ~ lim e _t x'(s) = ±c lim ( sec(x/2) + tan(x/2)) cos(x/2) = ±2c. 

Thus for each q e K. there exists a unique smooth solution v q of (3.2) on M + with 
v g(0+) = an d ^(0+) = noting that if q = then = identically. From (3.3), 
v£(0+) = for all q. 

Now let a PjQ denote the unit field with u = 9 + p and v = v q . The translational 
invariance o r p +2ir,q = & P ,q is refined by the glide symmetry cr p+7r g = cr p _ g , so the a p ^ q 
are in fact parametrized by an open Mobius band. In all cases, t 1— > t£,3 is an integral 
curve of cr p , q - Further qualitative aspects of the flow of a Pjq may be illustrated by the 
following two cases. We refer to the direction of £3 (resp. —£3) as "up" (resp. "down"). 

Case 1. p = 0. The sheaf of planes with axis L is invariant. Let a = \q\/2. If q > 
then the integral curves diverge from L inside the cylinder C a with equation r = a, 
moving upwards until cutting C a orthogonally, after which they move downwards and 
escape. The overall picture resembles a fountain. If q < the flow diagram is obtained 
by reflection in P and reversing the direction of flow. 

Case 2. p = tt/2. The family of coaxial cylinders C r (r > 0) is invariant. If q > 
then inside C a the streamlines are right-handed helices which spiral upwards. The pitch 
p r of the helices on C r diverges to infinity as r — * and converges to zero as r — > a—. 
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Outside C a the integral curves are downward-spiralling left-handed helices, with p r — > 
as r — > ct+ and p r — > oo as r — > oo. On C a the integral curves are horizontal anti- 
clockwise circles (when viewed from "above" ) . If q < the flow diagram is again obtained 
by reflection in P and reversing the flow, resulting in a family of upward-spiralling left- 
handed (resp. downward-spiralling right-handed) helices inside (resp. outside) C a . 

Regarding the bending of cr p , q , it follows from Proposition 1 that: 

i,-, i9 / / \9 sin 2 ^o(^) 
|Va M | 2 = v' q {r) 2 + \ 

which by a double application of l'Hopital's rule converges (to 2q 2 ) as r — > 0, in contrast 
to Examples 1 and 2. Since i>q(r) — > as r — > oo, \Va p ^ q \ 2 in fact remains bounded over 
IR 3 . However the total bending of a P)Q over M 3 is not finite. 

In all cases, when a PyQ is regarded as a map to S 2 its image is 5" 2 \ {—^3}, in contrast 
to Example 2 where the corresponding maps are all surjective. 

A more complicated pendulum equation, with damping and variable gravity, was used 
by R.T. Smith in [16] to construct "joins" of harmonic maps between spheres. 



4. New examples in hyperbolic space 

Now let M be a domain in H 3 . We denote the hyperbolic metric by ( , ), and initially 
work with the upper half space model, where ( , ) is: 

ds 2 = (dx 2 + dy 2 + dz 2 )/z 2 . 

We say that a unit vector field a is horospherical if o is tangent to a foliation T of M 
by horospheres. Without loss of generality, let T be the standard horosphere foliation, 
with leaves z = const. Let (£1, £ 2) £3) be the standard global orthonormal tangent frame: 

£1 = zd x , £ 2 = zd y , £ 3 = zd z , 

where d x = d/dx etc. Then (£1, £ 2 ) is a global orthonormal frame of J 7 , and the integral 
curves of £3 are the unit speed geodesies orthogonal to T. None of the £j are Killing. 
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The covariant derivatives of the are: 



£3, i = 1, f 6, i = 2, 

0, otherwise; t 0, otherwise; 



(4.1) 
-6, » = 1,2, 



0, i = 3. 

It follows from (4.1) by direct computation that the £j are harmonic, with constant 
bending: 

|Vei| 2 = l = |Ve 2 | 2 , |V6| 2 = 2. (4.2) 

The harmonicity of £j was observed in [9], in the general context of invariant unit fields, 
where it was also noted that none of the & are harmonic maps H s — > TH 3 . This may 
also be seen directly, by inspecting the additional equation [13] for a harmonic unit field 
a (on any manifold) to be a harmonic map: 

Y^R(<r,V Ei <T)Ei = 0, 

i 

where {Ei} is any local orthonormal tangent frame. For a non-Euclidean space form this 
yields the characterization that a harmonic unit field a is a harmonic map if and only if 
a is geodesic and solenoidal. Now £3 is geodesic but not solenoidal, whereas £1 and £2 
are solenoidal but not geodesic. In fact, since the horosphere foliation is not ruled, no 
horospherical unit field can be a harmonic map. 

Any horospherical unit field a has the following standard form: 

a = cos-u£i + sin-u£ 2 , (4-3) 

where u is a smooth locally-defined M-valued function. 

Proposition 2. Let a be a horospherical unit field, expressed in standard form (4.3). 
Then the bending of a is: 

|Va| 2 = 1 + |Vw| 2 , (4.4) 
and a is harmonic if and only if the following two equations hold: 

Au = 0, (4.5) 
u x sinu = u y cos u, (4-6) 
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where A is the hyperbolic Laplacian, and u x = du/dx etc. 
Remark. The hyperbolic Laplacian in the upper half space model is: 

-Au = z 2 (u xx + u yy + u zz ) - zu z , 

which is of course a linear second order elliptic differential operator. The nonlinearity of 
(1.1) appears in the first order constraint equation (4.6). 

Proof. It is convenient to introduce the following almost complex structure J in T: 

Jii = 6, Jii = -a- 

Then: 

V x cr = (X.u)Ja + cos-uV x £i + sin-uV x ^2, 
and it follows from (4.1) that: 

V^cr = (^i.u)Ja + cosw.^3, = {i"2-u)Ja + sinw^ 3 , ^£ 3 a = (£,3- u )J (T - 

Therefore: 

|Va| 2 = 1 + ^(Vw,&> 2 = 1 + \Vu\ 2 . 

i 

We now apply (2.1)-(2.3) to (4.3): 

V*Va = coswV*VCi + sinuV*Vf 2 

" 2 V Vcosu Ci - 2V Vsinu 6 + A(cosu)£i + A(sinu)£ 2 

= cosw^i + sin-u^ 2 + 2sin-uV Vu ^i - 2cos-uV v ^ 2 

— sin-uAw^i + coswAm^ + cosw|V-u| 2 £i + sin-u | Vw| 2 ^2, 

and then apply (4.1) and (4.4): 

= (l + |Vit| )<j + AuJa + 2(^i. u) sin u £3 — 2(^ 2 -u) cos-u^ 
= I Vcr | 2 cr + AuJa + 2z(u x sin u — u y cos u) £3. 

The result follows on comparison with (1.1). □ 
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One consequence of the constraint equation (4.6) is that, in contrast to the Euclidean 
case, if a is a harmonic horospherical unit field then the natural circle-action: 



e lt .a = cos(-u + + sin(u + t)£ 2 

generates a harmonic variation of a if and only if u is independent of x and y; other- 
wise, the only harmonic members of this S^-family of horospherical unit fields are ±cr. 
For example, by analogy with the Euclidean case (Example 1), let 9 be the following 
horospherical unit field, defined on M = H 3 \ L where L is the z-axis: 

^ = / O 9 ^ ^ / 9 9 

^x z +y z ^x z +y z 

Then 6 is a unit Killing field, and therefore harmonic [11]. However, apart from —9, the 
hyperbolic analogues of the other members of the 1-parameter family of Example 1 are 
not harmonic. In particular, the hyperbolic analogue of f is not harmonic. This does not 
contradict the results of [3] , since f is no longer radial. 

Theorem 1. Let a be a harmonic horospherical unit field, in standard form (4.3). 

(i) If u = u(x, y) then either a is invariant, or a = 6 upto translation in H 3 . 

(ii) If u = u(z) then a is a member of the following 2 -parameter family : 

o~p, q = cos(pz 2 + q)£i + sm(pz 2 + q)&, where p, q G K. 



Proof. 

(i) In this case (4.5) reduces to the Euclidean Laplace equation on R 2 . Let v.M 2 — > R 
be the harmonic conjugate of u. Then F = u + iv is a holomorphic function of ( = x + iy. 
By the Cauchy-Riemann equations: 

dF 

— = F C = U x ~ lUy, 

and hence: 

F^e tF = (u x - iu y ) e i ^ u+ ' lv "> = e~ v (u x - iu y )(cosu + i sin it) 
= e~ v (u x cosu + u y sinw + i(u x sin-u — u y cosu)). 
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Therefore (4.6) implies that the holomorphic function F^e lF is purely real, and hence 
constant. Consequently we obtain the first order holomorphic ODE: 

= ke~ lF , keR. 

d( 

This has solution: 

F(C) = -i\og(ik( + a), «GC, 

with real part: 

u(C) = &rg(ik( + a). 

If k = then u is constant and a is therefore invariant. Otherwise, by applying a 
parabolic isometry of H s , we may take a = 0, in which case a = 6. 
(ii) In this case equations (4.5) and (4.6) reduce to: 

and it follows that: 

u = pz 2 + q, p, q G R. □ 

The family a PjQ of Theorem 1 is 27r-periodic in q, and hence parametrized by a cylin- 
der; in fact, q is the parameter associated to the natural circle-action. This family of 
harmonic unit fields is interesting for a number of reasons. It follows from Proposition 2 
that amongst horospherical unit fields those which are invariant minimize energy (over 
relatively compact domains). These invariant fields are parallel when restricted to any 
horosphere, with respect to the induced (Euclidean) metric, and the same is true of the 
(Tp g. However, by Proposition 2, the bending of a PjQ is: 

|V(j Pi(? | 2 = 1 + |V-u| 2 = 1 + z 2 (u 2 x + u\ +u 2 z ) = l + 4p 2 z 4 . 

It follows that no member of the 1-parameter family a t = a t ,o is invariant under any 
discrete subgroup of hyperbolic isometries with compact quotient (except when t = 0) , 
but their local energies are arbitrarily close to the minimum. At first sight, the existence 
of the harmonic variation at appears to contradict the variational characterization of 
harmonic sections. However on a non-compact manifold harmonic sections are critical 
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points of the vertical energy functional with respect to all compactly supported smooth 
variations through sections, and the variation field of cr t is: 

— z 2 sin(tz 2 )£i + z 2 cos(t2 2 )^2, 

which clearly does not have compact support. 

5. Instability of H-parallel unit fields. 

In general, a harmonic unit field a is said to be unstable if there exists a compactly 
supported energy-decreasing variation of a. By [19] the Jacobi operator for a is: 

J a (a) = V*Va - |Va| 2 a - 2(Vff,V«)a, 

for all compactly supported smooth vector fields a pointwise orthogonal to a. The second 
variation or Hessian of the energy /bending functional at a is then: 

H„{a,P) = f (Ma),p)dV, 
Jm 

where dV denotes the Riemannian volume element, and a is unstable if and only if there 
exists a with H a (a, a) < 0. If U C M is a domain containing supp(a) then we say that 
a is unstable, or stretchable, in U . 

Theorem 2. Let a be an H-parallel unit field on H 3 . Then there exists R u > such 
that a is stretchable in any geodesic ball of radius R > R u . 

Proof. Without loss of generality we may assume that a = £3. Consider the unit field 
£ = £1, and let / be a smooth function on H 3 with compact support. Then a = /£ 
is a smooth compactly supported variation field for a, and since |Va| 2 = 2 by (4.2), it 
follows that: 

(Ma), a) = (V*Va,«) - 2/ 2 . 
Now V*V£ = £ and it follows from (2.1) that: 

( V *Va,a> = (V*V(/£),/£> = </V*V£ - 2V v/ £ + (A/)£, /£> = f 2 + f Af. 

Therefore: 

U»,a> = f Af - f 2 , 
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and hence by Stokes' Theorem: 

H a (a,a) = [ (\Vf\ 2 -f 2 )dV. (5.1) 

J M 

Now for any x$ G H 3 let Br(xq) be a geodesic ball about x$ of radius R. Fix 5 > 
and suppose that / is a bump function <p with: 



(p(x) = 



1, £ G Br(xo), 

0, xeM \B R+S (x ). 



Suppose further that <p is radially symmetric about xq, and piecewise linear with |V<p| = 
1/8 in the spherical shell Sr^r+^xq) = Br + s(xq)\Br(xq). Such <p may be approximated 
smoothly, with uniform convergence to both <p and V<p. For further computational 
simplicity we now work with the unit ball model of H 3 , where ( , ) is: 

2 _ 4(dx 2 + dy 2 + dz 2 ) 
ds ~ (1-r 2 ) 2 ' 

with r 2 = x 2 + y 2 + z 2 . We also take xq = (0, 0, 0). Then in terms of geodesic spherical 
polar coordinates (p, 6, <p) about xq we have: 

p = tanh(r/2), 

and: 

{1, O^p^R, 
R — p 

1 H — , i? < p < -R+ 1, 

0, p^R + 1. 
Let Vr (resp. Vr^+j) denote the volume of Br (resp. Sr : r + $) : and note that: 

<iV = sinh p sin dp dOdcf). 

It then follows from (5.1) that: 

n a {a,a) = ^ f dV ~ T2 [ (R + 5-p) 2 dV-f dV 
l-(R + S) 2 2(R + S) 1 

- ^2 ><R,fl+« - V R H £2 7l ~~ ^2 72 ' 
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where: 



7T r 2ir r R+S 



/ pdV =11 I psinh^ p sin <j)dpdOd<j) 

Js RiR+s JO JO JR 

"R+S 



2n / p(cosh2p — l)dp = 7r [psinh2p — icosh2p — p 2 ] 
.//,' 



21 i?+<5 
R 



and: 



h = p 2 dV = 2n p 2 (cosh2p-l)dp 

J Sr, r+s JR 

(p + |) sinh2p — pcosh2p 



= 7T 



We now note that the volume of a hyperbolic ball of radius p is: 



V p = 7rsinh2p — 27rp, 



from which: 



I x = (R + S)V R+S - RV R + ^cosh2 J R - J cosh(2i? + 26) + ttS(2R+S), 

7 2 = (i? + (5) 2 VR +5 - i? 2 VR + \ Vr, r+6 + ivRcosh2R - ir(R + 5) cosh(2i2 + 25) 
+ n5(4R 2 + 4RS + |5 2 + 1). 

Therefore: 

2(R + 5)h-I 2 = (R + 5) 2 V R+S - R(R + 25)V R - \ V R , R+5 

+ 7r<5cosh2i2 + nS(2RS + §5 2 - l), 



and so: 



Wa(a,a) = (Jj-^jv R , R+ s + ^cosh2i? + | (2R5 + §5 2 - l) . 



Finally: 



Vr,r +s = Vr+s-Vr = 7r(sinh(2i? + 2(5) -sinh2i?) - 2tv(R + 5) + 2nR 
= 27rcosh(2i? + 5)sinh5 - 2nd, 
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so that: 

HJa,a) = ^(2-8 2 )smh8cosh(2R + 8) + ^cosh2i? + % (2R8 + U 2 - 3) . (5.2) 



It follows from (5.2) that if 8 is sufficiently large then Ha (a, a) < for all R > 0. Let 8 U 
be the infimum over all such 5, and take R u = 8 U . □ 

Remarks. 

(1) From (5.1) we could have estimated: 

H a {a,a) ^ j \V V \ 2 dV - j ip 2 dV = ^V R , R+S - V R . 

However, in contrast to Euclidean geometry, the rapid rate of hyperbolic volume growth 
renders this upper bound positive for all R, 5 > 0, leaving us no alternative but to 
explicitly integrate (p 2 over the shell Sr jR+ s. 

(2) Further qualitative investigation of (5.2) shows that if 5 is sufficiently small then 
H a (a, a) > for all R > 0. Let 5 S be the supremum over all such 5. For all 5 G (8 S , 8 U ) 
there exists Ro > such that Ti a (a, a) < for all R > Rq. Numerical investigation 
yields estimates 8 S xs 1.471007 and 8 U xs 1.612195, and also shows that the transition to 
instability is extraordinarily rapid; for example if So = 1.471008 then R ^ 8.198206. 

(3) It is not known whether the energy functional on unit fields is minimized (over 
relatively compact domains) by the invariant horospherical fields, or even whether these 
fields are stable. 
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